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MEANS OF A DIRICHLET PROCESS AND MULTIPLE 
HYPERGEOMETRIC FUNCTIONS 1 

By Antonio Lijoi and Eugenio Regazzini 2 

Universitd degli Studi di Pavia 

The Lauricella theory of multiple hypergeometric functions is 
used to shed some light on certain distributional properties of the 
mean of a Dirichlet process. This approach leads to several results, 
which are illustrated here. Among these are a new and more direct 
procedure for determining the exact form of the distribution of the 
mean, a correspondence between the distribution of the mean and the 
parameter of a Dirichlet process, a characterization of the family of 
Cauchy distributions as the set of the fixed points of this correspon- 
dence, and an extension of the Markov-Krein identity. Moreover, an 
expression of the characteristic function of the mean of a Dirichlet 
process is obtained by resorting to an integral representation of a 
confluent form of the fourth Lauricella function. This expression is 
then employed to prove that the distribution of the mean of a Dirich- 
let process is symmetric if and only if the parameter of the process 
is symmetric, and to provide a new expression of the moment gener- 
ating function of the variance of a Dirichlet process. 

1. Introduction. The connections between properties of functionals of a 
Dirichlet process and the Lauricella multiple hypergeometric functions have 
been investigated in independent papers by Kerov and Tsilevich (1998) and 
Regazzini (1998); they also represent the point of departure of the present 
paper. The approach undertaken here is quite different from that of recent 
contributions to the study of exact distributions of functionals of a Dirichlet 
process. See, for example, Regazzini, Guglielmi and Di Nunno (2002). Some 
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specific properties of multiple hyper geometric functions are extended in such 
a way as to become significant properties of the laws of functionals of random 
measures. This is an unusual application of the theory of special functions. 
On the other hand, since these extensions can be thought of as infinite- 
dimensional versions of some fundamental types of special functions, the 
results presented in the following sections might be also of some interest 
from an analytic point of view. 

The main reason for analyzing the interplay between multiple hypergeo- 
metric functions and laws of Dirichlet functionals is twofold: it allows sim- 
plifications in the proofs of some well-known propositions and, more im- 
portantly, it leads to new results concerning the distribution of the above- 
mentioned functionals. Ties between the theory of multiple hypergeometric 
functions and some applied problems in statistics had been investigated in 
previous contributions, in addition to the papers already cited above. Dickey 
(1968) extends an identity due to Picard and applies it to Bayesian inference 
about multinomial cell probabilities, with a prior expressed by the Savage 
generalization of the Dirichlet distribution. Hill (1977) exploits suitable rep- 
resentations of Appell's hypergeometric functions to exactly evaluate the 
posterior moments of parameters of interest for the inference about vari- 
ance components. In a problem of Bayesian statistical inference for missing 
data, Dickey, Jiang and Kadane (1987) obtain representations of posterior 
moments and predictive probabilities in terms of ratios of Carlson's hyperge- 
ometric functions. Jiang, Kadane and Dickey (1992) explore computational 
methods for hypergeometric functions arising in Bayesian analysis. More- 
over, an introduction to Carlson's functions for statisticians can be found in 
Dickey (1983). 

As far as the present article is concerned, it is organized as follows. Some 
integral representations of the Lauricella hypergeometric functions, together 
with their probabilistic interpretations, are illustrated in Section 2. The Fei- 
gin and Tweedie existence condition and the Markov-Krein identity are then 
jointly restated in Theorem 1 in Section 3. The proof is based on some clas- 
sical results concerning Lauricella's hypergeometric functions. Despite its 
pure analytic nature, it is as simple as the proof of Theorem 9 in Tsilevich, 
Vershik and Yor (2000). See also Theorem 2 in Tsilevich, Vershik and Yor 
(2001). Moreover, we prove the existence of a one-to-one correspondence be- 
tween the set of nonnull measures, on M, with given finite total mass a > 
and the set of all probability distributions (p.d.) of linear functionals of a 
Dirichlet process with parameter having total mass a. Among recent papers 
focusing on such a bijection — sometimes referred to as the Markov-Krein 
correspondence — are those of Diaconis and Kemperman (1996), Kerov and 
Tsilevich (1998) and Tsilevich, Vershik and Yor (2000). In Section 4, ele- 
mentary properties of gamma processes lead to an expression for the p.d. of 
the mean of a gamma process in terms of the p.d. of the mean of a Dirichlet 
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process. Moreover, it is also seen that the p.d. of the mean of a gamma pro- 
cess is a generalized gamma convolution and is, hence, infinitely divisible. 
This fact suggests a simple proof of the absolute continuity, with respect to 
the Lebesgue measure on M, of the mean of a Dirichlet process (see Propo- 
sition 2 and Remark 1). In Section 5, an extension of the Markov-Krein 
correspondence is deduced from well-known integral representations of the 
fourth Lauricella hypergeometric function. Two applications of this exten- 
sion are considered in Section 6: the first determines the exact form of the 
p.d. of the mean of a Dirichlet process, and the second is a characterization 
of the Cauchy distribution. Section 7 is devoted to a representation of the 
characteristic function of the mean of a Dirichlet process, via a multidimen- 
sional extension of Rummer's confluent hypergeometric function. The latter 
can be considered as a confluent form of the fourth Lauricella function and 
admits a representation as a contour integral proved by Erdelyi (1937). The 
distribution of a vector of means of a single Dirichlet process is examined 
in Section 8. The identities given in Section 5 are trivially extended to this 
setup. This is also helpful in proving the absolute continuity, with respect 
to the Lebesgue measure on M. d , of a vector of means of a single Dirichlet 
process, a finding which answers a question raised in Regazzini, Guglielmi 
and Di Nunno (2002). Finally, a representation of the moment generating 
function of the variance of the Dirichlet process is provided and it would 
seem a meaningful improvement on previous contributions to the subject. 
See, for example, Cifarelli and Melilli (2000). 

2. Probabilistic interpretation of the Lauricella fourth function Fn- The 

topic of multiple hypergeometric functions was first approached, in a sys- 
tematic way, by Giuseppe Lauricella at the end of the 19th century. See, for 
example, Exton (1976). He defined four functions which are named after him 
and have both multiple series and integral representations. In particular, the 
fourth of these functions, denoted by Fjj, admits integral representations of 
importance in our treatment. Given any £ = (£i,...,£ n ) in M n , set |£| for 
J2k=l £k an d (•, •) for inner product. Moreover, let T n : = {u = (u±, . . . , u n ) £ 
W 1 :U{> for every i and |u| < 1}. With this notation, an integral represen- 
tation of Euler type [see (26) in Lauricella (1893)] is 

r(6 1 )---r(b n )r(a-|b|) E1 

— F D (c, h, . . .,b n ',a;xi,.. .,x n ) 

r(o) 

(2.1) 

= / <- 1 ---<»- 1 (l-|u|) a -| b |- 1 (l-(u,x))- c ^ 1 .-- dll n , 

which is valid for every x = (x\, . . . , x n ) in [0, l) n and h = (b\, . . . , b n ) with 
strictly positive real part, that is, Re(ftfc) > for every k, and such that 
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Re(a — |b|) > 0. A further representation provided by formula (25) in Lau- 
ricella (1893) is 



Fd(c, bi,...,b n ;a;xi,...,x n ) 



(2.2) 



[0,1] 



[1 - UXl) 



(1 — ux n ) b " B ( du; c, a — c) 



which holds true whenever Re(c) > and Re (a — c) > 0, with 

(5i(A), ifRe(a-c) = 0, 

r(a) 



B(A;c,a- c) = < 



F(c)T(a-c 

x 



u^ 1 (1 - u)^" 1 du, if Re(o - c) > 0, 

'An(o.i) 

for every A in the Borel class on R, 3S(M), and 6 X standing for the unit point 
mass concentrated at x. Hence, letting a = c, one gets 

(2.3) Fn(a,bi, . . .,b n ;a;xi, ...,x n ) = (l- x\)~ bl • • • (1 - x n )~ bn . 

Representation (2.1) has an obvious probabilistic interpretation, which was 
first stressed and exploited by Carlson (1977), giving rise to what he called 
"method of Dirichlet averages," for unifying a part of the theory of special 
functions. In point of fact, let <p be a random probability measure supported 
by S = {0,xi, . . . , x n }, with (xi, . . . , x n ) in (0, l) n , and assume that the ran- 
dom vector (<p{x\}, . . . , <p{x n }) has the Dirichlet distribution, @, with pa- 
rameter (6i, . . . , b n , a — |b|), and a — |b| > 0, b^ >0 for every k, that is, the 
distribution characterized by 

T(a) 



9{d-a) 



r(6 1 )---r(6 n )r( a -|b| 



X U 



bl-1 



■ib-Hi 



u 



sa-|b|-l 



It„(u) du± ■•■ du n , 



where Ib denotes the indicator function of set B. Whence, combination of 
(2.1) with (2.3) gives 



(2.4) 



1 



x(p{dx) 



exp 



log(l — x)a(dx) >, 



provided that a is the measure on the power set of S determined by a{xk} = bk 
for k = 1, . . . , n and a{0} = a — |b| . Here, and in the sequel, given any proba- 
bility measure tp and real- valued function h such that f \h\ dtp < +oo, <p(h) 
denotes J hdp. 

Equality (2.4) represents the most elementary version of an identity estab- 
lished by Cifarelli and Regazzini (1979a, b, 1990) for functional Dirichlet pro- 
cesses with parameter a. Recall that, given a finite measure a on (M.,£3(M)) 



DIRICHLET PROCESS AND HYPERGEOMETRIC FUNCTIONS 



5 



with a := a(R) > 0, a random probability measure (p is said to be a (func- 
tional) Dirichlet process with parameter a if, for every finite and measurable 
partition C±, . . . , C n , C n+ i of R, the random vector (tp(C\), . . . , (p{C n )) has 
the Dirichlet distribution with parameter (a(Ci), . . . ,a(C n ), a(C n+ i)). As 
for the definition and the main properties of a functional Dirichlet process, 
the seminal contribution in Ferguson (1973) still represents a sound refer- 
ence. 

Throughout the following sections, C( w \ a > f) wm indicate the integral 
/ K log(l + wf(x))a(dx), while logz will denote the principal determination 
of the logarithm of the complex number z, that is, logz = log |-z| + i arg(z), 
where arg(z) is chosen in (— tt, tt]. 

3. Lauricella theory and Markov Krein identity. Diaconis and Kemperman 
observed that the identity of Cifarelli and Regazzini cited above is closely re- 
lated to a well-known version of the Markov moment problem. See Akhiezer 
and Krein (1962), Diaconis and Kemperman (1996) and Kerov (1998). As 
a matter of fact, some authors refer to this identity as the Markov-Krein 
identity. Compare, for example, Kerov and Tsilevich (1998) and the more 
recent paper by Tsilevich, Vershik and Yor (2000). According to Kerov and 
Tsilevich (1998) and Regazzini (1998), the Markov-Krein identity extends 
(2.4) from the (finite-dimensional) Dirichlet distribution to the functional 
Dirichlet distribution. Notice that in the real domain the extension holds 
true provided that the support of the parameter of the latter distribution 
is bounded above. In the complex domain, a variant of the same identity, 
which holds for any parameter, is established in Kerov and Tsilevich (1998) 
and in Regazzini, Gugliehni and Di Nunno (2002). In the following section, 
the Lauricella theory is exploited to state jointly the Markov-Krein iden- 
tity and the Feigin and Tweedie condition for finiteness of the mean of a 
Dirichlet process. 

3.1. The Markov-Krein identity. Given any finite measure a on (R, &(M)) 
such that a(R) = a > 0, let 3> a denote the functional Dirichlet p.d. with pa- 
rameter a. A suitable sequence of simple functions £ n , n > 1, exists that 
satisfies £ n (x) T I{x) if x > 0, £, n (x) j I(x) when x < and |£ n | | |/|, where 
/ is the identity map on R. Therefore the sequence of finite measures a n 
defined by 

a n (B)=a{£ n eB}, Bef(I) 

converges weakly to a, namely a n =>■ a. Now define P to be the space of all 
probability measures on (R,^(R)) endowed with the topology of weak con- 
vergence, and let & be the Borel u-algebra on P. In this way, it follows that 
the identity map dp on (P, Si a ) is a random probability measure with p.d. 
3i a . Moreover, (p n := (p o has p.d. & an . A straightforward application of 



6 A. LIJOI AND E. REGAZZINI 

the Lauricella formula (2.4), combined with a standard analytic continuation 
argument, yields 

(3.1) / (l + it<p(foZ n ))- a $> a {d<p)=exp{-((it;a n ,f)}. 



Theorem 1. Setting L:={<pEF:cp(\I\) :=lim<p(|£ n |) is finite} , one has 

® a (L){l-9 a (L)} = U, 
and the following two conditions are equivalent: 

(i) 9 a {L) = l. 

(ii) Jlog(l + \x\)a{dx) < +oo. 
Moreover, ifS# a (L) = l, then 

L (i + ^(j))^ q( dip) = ex ^~^ J »' te M - 

The equivalence of (i) and (ii) was first proved by Feigin and Tweedie 
(1989), but it was already contained — at least in part — in Cifarelli and 
Regazzini (1979a, b). A complete development of the argument used in the 
latter paper can be found in Cifarelli and Regazzini (1996). Here it is shown 
that both equivalence of (i) and (ii) and identity (iii) easily follow from (3.1). 

Proof of Theorem 1. In view of the definition of (£ n )n>i, the oper- 
ations of integration and of taking limit can be interchanged to obtain 

((it;a n ,\I\) 

I" log(l + t 2 £(x))a(dx) + i [ arg(t|f n (z)|)a(dz) 



-i /loe(l + tV)a(<fc)-H / «rg(tM)a(*0, asn-+oo, 
Jr Jr 

where an infinite real part in the above limiting expression is allowed. More- 



over, 



[i + iM\UT a ^a(dip) 

= ^exp|-^log(l + iV(l£n|)) - iaarg(ty>(|£ n |)) W«( dtp) 



1 



L ( i + *M|/|))^ W ' as ^+°°' 

since = +oo if (p G L c . Hence, by (3.1) with /(•) = | ■ |, one gets 

(3.2) / (l + i^(|/|))- a ^(^)=exp{-C(it;a,|/|)} 
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for every t 6 R, with the proviso that the right-hand side is whenever the 
real part of Jlog(l + it\x\)a(dx) is infinite. Thus, if J R log(l + \x\)a{dx) = 
+oo, the right-hand side of (3.2) is for every t ^ and, therefore, @ a {L) 
must be 0. On the other hand, if / R log(l + \x\)a{dx) < +oo, then t i— ► 
exp{— J R log(l + it\x\)a(dx)} is continuous at t = 0, that is, 



limexp< — / log(l + it\x\)a(dx) > = 1, 



and, by taking the limit (as t — > 0) on both sides of (3.2), one obtains 

S a {L) = l. 

Finally, if (ii) holds, interchanging the operations of integration and of taking 
limit in (3.1), with f = I, yields (hi). □ 

Identity (hi) in Theorem 1 can be invoked to prove that the p.d. of (p 
characterizes the parameter of 3) a within the class of all measures a on 
(R,=^(R)) for which a(R) has some fixed value a > 0. This is the essence 
of the Markov-Krein correspondence and it is dealt with in Section 3.2. In 
addition to the papers listed in Section 1, Andrea Ongaro, in a personal 
communication, has shown us a proof of such a correspondence based on a 
completely different approach. 

3.2. Uniqueness theorem. Define F to be the set of all finite measures on 
(R,^(R)), and put 

F a = ja G F:a(R) = a and J log(l + \x\)a{dx) < +oo|, 

M a = {& a o{p(I)- 1 :ae¥ a } 

for every a > 0. In words, M a is the set of all p.d.'s of (f(I) when (p is a 
Dirichlet process with parameter varying in ¥ a . Let A stand for the distri- 
bution function associated to a, and A for (a — -A)I[o,+oo) — ^(-oo,o) • Clearly, 
to any a in F a there corresponds a unique fi a := S> a o (p{T)~ ■ Moreover, one 
has 

Theorem 2. Any a in ¥ a is determined by its fi a . 

Proof. Fix fi ai in M a and suppose 02 £ F a is such that fi a2 = fi ai . 
Then 

exp{-C(^; ai, I)} = / m 7J^-y a ^ (*«0 

(3 3) f 1 

= / n 1 \ a f la ^ dx ) [by hypothesis] 
Jr (1 + itx) a 

= exp{— £(it; ct2, 1)} [from (iii) in Theorem 1]. 
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Use integration by parts to obtain 

C(it;otj,I)= f * Aj(x)dx, teR,j = l,2, 
Jr 1 + itx 

which, combined with (3.3) and analytic continuation, yields 

(3.4) / ilM dx = ( dx, zeC with lm(z) / 0. 

Ju z + x Ju z + x 

Taking, l(z, ctj) for J tL [Aj(x)/(z + x)] dx, j = 1,2, and resorting to the Stieltjes- 
Perron inversion formula [see Theorem 12.10d in Henrici (1991)], one has 

= ^— lim{^(-C - ie, «i) - K~£ + »e, «i)} 
= — lim{Z(-£ - ie, q 2 ) - + ie, a 2 )} 

ZTTl ej.0 

provided that ^ is a continuity point for both A\ and ^2- This suffices to 
conclude that a,\ = «2, since Oi(R) = a, i = 1, 2. □ 

Theorem 2 states that there is a bijection /3 a of F a to M a . Clearly, M a C P 
for every a > and, more precisely, M a C P for each a > 0. Namely, for every 
a > 0, there is some probability \i on 3§(M) which is not the p.d. of <f(I) for 
any 3> a with a £ ¥ a . For example, take 

fJ>( dx ) = ^"{^-I-tj-I)^) + I (l,l+r?)( ;c )} 

with 77 > 0. It is easily seen that such a /x cannot be the p.d. of the mean of 
a Dirichlet process with parameter a, for any choice of a in F a and for every 
a > 0. One can argue by noting that the support of the p.d. of the mean of 
a Dirichlet process with parameter a must coincide with the closure of the 
convex hull of the support of a. This can be proved by combining a result 
concerning the topological support of the Dirichlet prior given in Majumdar 
(1992) and the equivalence of (i) and (ii) in Theorem 1. We finally remark 
that general properties of M a are discussed, for example, in Section 2.3 of 
Kerov (1998). 



4. Markov Krein identity and means of gamma processes. From ele- 
mentary properties of the gamma process it follows that the right-hand side 
of (hi) in Theorem 1 is the conjugate of the characteristic function of the 
"mean" of a gamma process with parameter a. See (4.1). It is easy to ver- 
ify that such a characteristic function coincides with the Fourier-Stieltjes 
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transform of a generalized gamma convolution, according to the terminol- 
ogy introduced by Olof Thorin. See, for example, the original contributions 
in Thorin (1977a, b, 1978a, b) and the systematic treatment in Bondes- 
son (1992). From one of the Thorin results, it is possible to state the infinite 
divisibility of the distribution of the mean of a gamma process. This fact 
motivates the search for the Levy-Khintchine representation of the above 
characteristic function. And this representation, in turn, allows us to deduce 
immediately the absolute continuity, with respect to the Lebesgue measure, 
of the corresponding distribution. 

4.1. Remarks on the characteristic function of the mean of a gamma pro- 
cess. Endow F with the topology of weak convergence, and let a be an 
element of F a , with a > 0. Denote the Borel a-field on F by and call 
functional gamma distribution with parameter a the probability measure 
^ a on (F,^") defined as follows. Set 7 for the identity map on F and say 
that 7 has the functional gamma distribution with parameter a, @ a , if for 
every finite and measurable partition {B\, . . . of R, the random vari- 
ables 7(-Bi), • • • , 7(-Bfe) are independent, with gamma distribution such that 
-E(7(-Bj)) = Var(7(5j)) = a(Bj), for each j = 1, . . . , k. It is well known that 
7(-)/7(R) is a random probability measure with p.d. S> a . Moreover, one 
easily obtains the following representation for the characteristic function of 
7(1) =fayy(dx): 



Thus, according to Section 3.1 in Bondesson (1992), the p.d. of the random 
mean is an extended form of generalized gamma convolution, with 

Thorin measure a*, where a* = aoj^ 1 and J(x) = l/x for any x in R\{0}. 

On the other hand, (4.1) can be invoked to establish a representation of 
the p.d. of 7(1) in terms of the p.d. fx a of (p{I). 

Proposition 1. Suppose that a is an element in ¥ a , for some a > 0, 
and < a{0} < a. Then, the p.d. 0/7(7) is absolutely continuous with respect 
to the Lebesgue measure and 



(4.1) 



teR. 





is a density function of such a p.d. 
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PROOF. Theorem l(iii) yields 




z a-i e -z(l-itx) dz ^ dx ) 




r(a) J (-00,0) 



x 



a e -y/x 



li a {dx) dy 



+ 




X 



a e- y/x n a {dx) dy 



where the last equality follows from the application of the change-of- variable 
formula, with y = y(z) = zx, and of the Fubini theorem on iterated integrals. 
It follows that exp{— £(— it; a, I)} =^ a (e^W) is the Fourier transform of 
the probability density function q(-). □ 

Absolute continuity extends to p a . We propose a new and simpler proof of 
this statement. See Regazzini, Guglielmi and Di Nunno (2002) for a different 
line of reasoning. 

Proposition 2. Suppose that a is in F a , for some a > 0, with < 
a{x} < a for every xel. Then, the p.d. of <p(T) under 3> a is absolutely 
continuous with respect to the Lebesgue measure on R. 

PROOF. Suppose that fj, a has a singular part ^ a>s such that fx a ,s[c, c + 
A] > with A > and c in R. Since a is nondegenerate, the p.d. function 
G c a of 7(7 — c) is absolutely continuous and, by virtue of Proposition 1, 



G£(A) - G£(0) = (r(a)r 1 / A x a ~ l / (0)+oo) y-'e-B/VfiZidy) dx, fx% being the 



p.d. of ip(I — c). Without real loss of generality, it suffices to consider the 
case c = 0. Then, for any finite class of disjoint intervals (pi,p'i), • • • , (Pm,Pm) 
with px > 0, P2 > p'i, ■ ■ ■ , Pm > Pm-i ' an d f° r an y 5 > 0, one has 



£[G a ((l + <^)-G Q (p,-)] 
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rl+S 1 m 

Because of the existence of the singular part p a)S (with c = 0), there is e > 
such that: With each r/ > one can associate a finite class of disjoint in- 
tervals, included in [0, A], (pi, p[), . . . , (p m , p' m ) such that J2(p'j ~ Pj) < V 
and, for some e > 0, Pa[Pj, Pj] > £• On the other hand, from the abso- 
lute continuity of G a , there is r] such that J2{G a ((l + — G a (pj)} < 

any finite class of disjoint intervals satisfying 
J2(p'j ~ Pj) < ?7> yielding a contradiction and, thus, completing the proof. □ 

4.2. Levy-Khintchine representation o/exp{— £}. The next proposition 
extends part of Theorem 3.1.1 in Bondesson (1992). It involves the distri- 
bution function A associated with a, the function g defined by 



(0,+oo) 

x E 



Ho,+oo)(x) I e yx dA(y 1 

1 ~r X L J(0,+oo) 

(4.2) 

+ I(-oo 1 o)W / e-^dAiy- 1 

J (-oo,0) 

and the well-known fact that the Levy-Khintchine representation of an in- 
finitely divisible characteristic function, determined by the pair (7, G), is 
expressed by 

i_ exp /< 7 t+ / ( e itu_ 1 _Jp_\tt^ dG{u 



1 + u 2 J u 

where 7 S R, and G is nondecreasing, right-continuous with G{— 00) = and 
G?(+oo) < +00. 

Theorem 3. Suppose a satisfies condition (ii) in Theorem 1. TTien £/ie 
characteristic function of 7(1) is infinitely divisible with Levy-Khintchine 
representation determined by the pair (7,6?), where 7 = J K \{ } M" 1 ^^) ^ x 
and dG(x) = g(x) dx for every x G M, g freing defined by (4.2). /n particular, 
^(e^W) = 1 if and only if a = 5q. 

Proof. The result is first proved for the case / R |a;|a(dx) < +cxd. By 
differentiation under the integral sign one gets 

d , , i«(rv /" ix 



4log^(e^)= / -^-a(rfx) 



+00 



is / e- 2(1 -* te) ^^(s) 
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-i [ °° [ e-^dAiy- 1 )) d£ 

Jo W(0,+t») / 

+ 1 [ f e-^dA(y- 1 )) d£, 

J-00 \J(-oo,0) / 



'(-oo,0) 

where the application of the Fubini-Tonelli theorem is valid in view of the 
above extra-assumption. Next, by defining g as in (4.2), one obtains g > 
0, J R g(x) dx < +00, / K \x\~ 1 g(x) dx < +00. Hence, G{x) := J_f g(u) du is 
well defined for every x in R, and it turns out to be nondecreasing with 
G(— 00) = and G(+oo) < +00. Moreover, letting 



4>*(t):= l_(e itx -l) 1 -^g(x)dx, t€ 



then e^* is an infinitely divisible characteristic function with 

7 = / \x\~ 1 g(x) dx 



and dG{x) = g(x) dx, such that 

^(t) = ^iog%( e ^), teR. 

The latter can be verified by interchanging the derivative with the integral 
in the expression of </>*(•)• It implies that is the characteristic function 
of 7(1) when J R |x|a((ix) < +00. 

The extension of this conclusion to any a such that J R log(l + \x\)a{dx) < 
+00 can be carried out through direct calculation of (/>*. □ 



Remark 1 . A nice application of this theorem is a further proof of the 
absolute continuity of the p.d. of In fact, under the hypotheses in 

Proposition 2, one has Jr\{ } x~ 2 g{x) dx = +00 and this, applying a crite- 
rion due to Tucker (1962), suffices to state the absolute continuity of the 
(infinitely divisible) p.d. of 7(1). 

5. Extension of the Markov Krein identity: the Lauricella identity. In 

Section 3, the so-called Markov-Krein identity has been obtained as an 
extension of the Lauricella formulae (2.1) and (2.2) with a = c > 0. The 
Lauricella general representations are now employed to prove an identity 
that holds for all pairs (a,c) of strictly positive numbers. For the remainder 
of this section, 

/•(!+) 

1\ := / exp{— ((iwt; a, I)}b(w; c, a — c) dw 
Jo 
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is to be meant as the integral along the contour OGG' , as shown in Figure 1, 
of the function exp{— £(iwt; a, I)}b(w; c, a — c), where 

r» 



b(w; c,a — c) 



w) 



a— c— 1 



T(c)T(a-c) 

for every w in C such that Re(u>) > and Im(u>) ^ if Re(iy) 6 (0, 1]. 

In view of the Cauchy integral theorem, this contour can be deformed into 
the path of integration consisting of: (A) the straight line segment w = x — ir 
with x varying in (0, 1 — e), (C) the circle 1 + ee ld with G {—it + r/,7T — rj), 
for a suitable rj > and (A') the straight line segment w = x + it. Thus, if 
c and (a — c) are strictly positive, with (a — c) different from 1,2,..., one 
can let e — ► to obtain 

(5.1) ii = (1 - e 2 ^ (a ~ c) ) / exp{-C(i^; a, /)}&(«;; c, a -c)dw. 



This is required to state the following extension of point (hi) in Theorem 1, 
which we shall call Lauricella' ' s identity. The symbols B and fj, a appearing 
in this identity have the same meaning as in Sections 2 and 3, respectively. 

Theorem 4. If, for some a > 0, a is in ¥ a , then 
(5.2) / (1 + itx)~ c ii a (dx) = / exp{— ((iut; a, I)}B( du; c, a — c) 

JR J [0,1] 

if a>c> and 

(1 + itx)~ c fj, a (dx) 

r( c -a + i)r(o) 



(5.3) 

ifc>a>0 



2vrir(c) 



exp{— ((iwt; a, I)}w c (w — 1 



a— c— 1 




Fig. 1. 
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PROOF. First note that, if < c < a, the following equality holds true 
by virtue of analytic continuation of both sides of (2.1) and (2.2): 

/ t-\ i u~ it = / exp{-({iut;a n ,I)}B(du;c,a-c). 

Jp (1 + ity(£ n )) c 7[0,1] 

Next, argue as in the proof of Theorem 1 and apply the basic limit theorems 
of integration theory to obtain (5.2). 

Finally, suppose that a > c > and use what has just been proved and 
(5.1) to state the equality 



/ 



(l + itx) c ^ ay ' 1 — e 2Tvi(a-c) 



when a — 1,2, ... . Now, take Euler's reflection formula r(^)r(l — z) 
7r(sin7rz) _1 to obtain 



:Ha(dx) 



(1 + itx) c 

r(c-a + i)r(o) . i, 

w c - l {w - if - 1 exp{-((iwt; a, I)} dw. 



2vrir(c) 



Notice that the right-hand side of the previous equality, as a function of 
c, is analytic on {c 6 C : Re(c) > 0, Re(a — c) ^ {0, 1,2,.. .}}. Hence, such an 
equality can be extended to any c with c > a. This completes the proof. □ 



The Lauricella identity establishes that the Stieltjes transform (of order 
c > 0) of fi a can be viewed as a mixture of the Stieltjes transforms (of 
order a) of p.d.'s of means of Dirichlet processes with parameter a w (-) := 
a{x :wx G •}. Apropos of this, recall that 

r(c-a + l)r(a) 



/ ™ c -V-l) a dw = l. 
Jo 



Compare, for example, 3.1.27 in Slater (1960). 



6. Applications of the previous results. 



6.1. Exact forms of the p.d. of f x(p{dx). Much work has been done 
on the exact form of the p.d. of <f(I). As the following arguments show, 
Lauricella's identity considerably simplifies the solution of the problem. 
First, rewrite Theorem 4 with c = 1 and a such that a{x} < a(R) for 
every x. Denoting a density function for \x a by m a and setting L(w) := 
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Ir( w + x)~ 1 m a (x) dx, one has 

w~ l I ex.p{—((uw~ 1 ;a,I)}B(du;l,a — l), if a > 1, 
J[0,1) 

w -1 exp{— £(uj _1 ; a, I)}, ifa = l, 

^ _1 (l-exp{2vri(a-l)})" 1 
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L(w) = < 



(1- 



exp{— ((uw ;a,I)}b(u; l,a — 1) du, if a E (0,1), 



for every w E C for which Im(w) 0. Thus, the Stieltjes-Perron inversion 
formula yields 



U a (xi,X2\ = — limlm 

7T e|0 



■>:> 



1 



if a = 1, 

/ u a (xi,x 2 ] 

a- 1 



lim Im 



7T e|0 



Xl (-A - ie) 



1 



exp<^ -C 



1 



-A — ie 



a, I > d\ 



Xl (-A - ie) 
exp. 

o 



-A — ie 



a,/) j(l-u) a " 2 du dX 



if a > 1 , and 

Ma(xi,X 2 ] 



— lim Im 

7T e|0 



X2 



Sl (1- e 2 ^( a - 1 ))(-A-ie) 



y exp|— — —;a,Ijjb(w;l,a — l)dwd\ 



-A — ie 

if a 6 (0, 1). In particular, when a = 1, one obtains 



m Q (£) 



1 



lim Im 



1 



7T e|0 (— £ - ie) 



-£ — ie' 



;a,I 



if a > 1 and the saltus of A at each discontinuity point is strictly smaller 
than 1, one has 



1 1 
m a (C) = - lim Im— — 

vr e|0 (-£ - «e) 7[o,i] 



exp<^ -C 



-£ — 



; a, I^j |.B( rfn; 1, a — 1). 



At this stage we can achieve the same results as in Proposition 9(h) and (iii) 
of Regazzini, Giglielmi and Di Nunno (2002), by a simple application of the 
change-of-variable formula. 
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Example 1. (a) Let a(dx) = Cg a i{ dx) be the Cauchy density function 
with parameters 9sR and a > 0, that is, 



a 



7r(l + a 2 (x-0) 2 ) 



dx. 



x 6 



It is also convenient to identify Cq :+oq with the degenerate distribution 8g. 
By resorting to the expression of m a given above when a = 1, it is easily 
verified that, for any 2; in R, 

m Q (x) = — exp{ — [ log \p — x\Cg a ( dp)\ sin( — h arctan(cr(a; — 6)) 

7T I JWL\{x} ' J \2 

a 



n(l + a 2 (x-e) 2 ) 



if £ (0 + 00), whereas p a is degenerate at if a = +00. In other terms, 
if a coincides with Cq a 2 for some 9 £M. and a £ (0, +00], then \i a = a, a 
well-known result stated by Yamato (1984). In Section 6.2 the converse will 
be proved. 

(b) If a(dx) = al^i) (a?) dx with a > 1, then, for any £ £ (0, 1), 



(a - 1)(1 -£) a_1 e a ^ 



7T 



e o(l-f)u u o-2 S i n ( a7r (i _ £)(i _ „)) 

[(1 - £)U + f]aM-a(l-£)u[(l - f )(1 - u )]a(l-S)(l-«) 



(c) If a(cfe) = a(cr\/27r ) 1 exp{- -^(x - 9) 2 } dx, for any a > 0, 6 £ R and 
a > 1, then application of 2.6.22.1 in Prudnikov, Brychkov and Marichev 
(1986) yields 



m a (£) 



(a -IX 



a-l 



7T 



(1-U) 



x exp 



a-2 



sm 



a 



owr^ 1 - $ 



£-u<9 



ctV2tt 
5 



cxp 



2u 2 cr 2 



z> 1 (£-m0) 5 



2'2 



2a 2 



u=l 



du, 



for any £ > 0, where <3? denotes the distribution function of a Gaussian ran- 
dom variable with zero mean and variance equal to 1 and \& is the Tricomi 
confluent hyper geometric function. On the other hand, letting £ < 0, one has 
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7T 



x / (l-u) a " 2 sin 



o 



air® 



x exp < -= exp 



2 



a f (f - «0) 



01/ V V2'2' 2a 2 



</=i 



6.2. Characterization of the Cauchy distribution. We prove the charac- 
terization of the Cauchy distribution already mentioned in Example 1(a). 
The proof is based essentially on the result in Theorem 4. As in Theorem 2, 
[3 a is the bijection of M a on F a . 

Theorem 5. The class of all fixed points of the bijection (3\ is {C$ a 2 :9 G 
M.,0- 2 G (0,+cx>]}. 

In other terms, if a G Fx, then fi a = a if and only if a is Cauchy or degen- 
erate. In Section 4 of Cifarelli and Regazzini (1990) another characterization 
is given under stronger hypotheses, that is: A nondegenerate probability is 
a fixed point of [3 a for every a if and only it is Cauchy. 

Proof of Theorem 5. In the light of Example 1(a), it suffices to prove 
that the condition is necessary. Suppose that fj, a = a G Fx. Then ii na = a for 
every n G N. Indeed, by the Markov-Krein identity, 

(6.1) / 1 . a(dy) = exp(- / log(s + iy)a(dy)\ 

Jm.s + iy I it J 

holds for every s£M \ {0} by virtue of the above hypothesis. Differentiate 
both sides to obtain 

7 — -^~a(dy) = exp(- f log(s + iy)a{dy) X [ ——a(dy) 
(s + iyY I Ji J Jr s + ty 

= exp|-2^1og(s + iy)a(dy)^ [by (6.1)]. 

Hence, the result holds true for n = 2. Assume that it is valid for all m <n. 
Then 



/ 7 r~zr a (dy) = exp( — (n — 1) / log 

Jr {s + iy) n 1 I Jr 



(s + iy)a(dy) 
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and, by differentiation, 

/ ( T~- \n a ( d y"> =exp(-(n- 1) [ log(s + iy)a(dy)\ [ 1 . a(dy) 

= expj-ra^log(s + iy)a(dy)\ [by (6.1)]. 

By induction one can conclude that fi na = a for every n > 1. Next, this 
statement and Theorem 4, with a = 2 and c = 1, give 

1 /" 1 Z" 1 

-fjL na (dx)= / exp{— £(zui; no, J)}(1 — u) n ~ 2 du 



n — 1 jm 1 + ifcc 
where, by hypothesis, 



/ 1 -TT Vna(dx) = [ /j, a (dx) = exp{-((it;a,I)}. 



Hence, 

exr> -f — t ( it', cv . TW = / exn-f — C( in: nn. TW\ 1 

n — 1 



— - exp{— ((it; a, I)} = J exp{— ((iu; na, I)} ( 1 — — J du. 
In particular, letting ((it; a, I) = p(t) and n = 2, one has 

te -p(t) = /" e -2p(«) dXf t G R _ 



Differentiation of both sides gives e p ^ — tp'(t)e p ^ = e 2 p^\ that is, 

p(t) = log(l + wt) 

for some w G C. Hence, exp{— J R log(s + ix)a(dx)} = (s + ui) -1 holds true 
for every s > 0, and by (6.1) one gets 

(5 + ^)-!= / (s + iy)- l a(dy)= j +C ° e~ sx ( f e-* v a(dyj) d£. 
Jr Jo \Jr J 

The left-hand side is the moment generating function of e~ wx l(Q^ +00 ^(x) for 
every s > —Ke(w), implying e~ wx = / R e~ lxy a(dy), for every x > 0. Analo- 
gously, e™ 1 = / K e~ txy a(dy) is valid for any x < 0. Therefore, 

a(dy) =C e>a 2(dy), j/GR, 

where a = l/Im(to), Im(w) > and 9 = ~Re(w). □ 

Remark 2. Richard Olshen has drawn the attention of the authors to 
the following problem. Let (X n ) n >\ be a sequence of exchangeable real- 
valued random variables, whose p.d. P is supposed to be characterized by 

(6.2) P(A 1 x • ■ • x A n x R°°) = / p(A 1 ) • • -p(A n )9 a ( dp), 
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n £ N, Ak S 3S(M) for k = 1, . .. ,n, a being any element of Fi. Then, by 
de Finetti's representation theorem and the strong law of large numbers for 
i.i.d. random variables, one has 



P 



Y,X k /n-(p{I) 



k=l 







<p 



which entails 



n 



k=l 



x(p{dx) 



1 a.s.-P, 



a.s.-P. 



Thus, the limiting p.d. of Y^k=i x k/ n coincides with the law of <p(I), that 
is, with fi a . Moreover, recall that each has p.d. a. A natural question 
is whether the above limiting p.d. may coincide with a. An answer can be 
deduced from Theorem 5, that is: If (X n ) n >x is a sequence of exchangeable 
random variables with a law characterized by (6.2), then the limiting distri- 
bution of the empirical mean J2k=i x k/ n is a if and only if a is an element 
of {C „ 2 :eeR,o- 2 e(O,+oo]}. ' 



7. Characteristic function of f x(f(dx). Let b\, . . . , b^ be strictly positive 
numbers such that |b| < a, and let x\,...,x n be arbitrary real numbers. 
Erdelyi (1937) defines a confluent form of the fourth Lauricella function, 
through the following limiting process: 

n <&(bi,...,b n ;a;itxi,...,itx n ) = ]haF D (e~ 1 ; h, . . . , b n ; a; ietx±, . . . , ietx n ), 

e|0 



where t is any real number. From 
of 3> a , we get 



.5) in Erdelyi (1937) and the definition 



r(a) 



r(6 1 )---r(6 n )r(a-|b| 



u 



bi-i 



■•^(l-lul) 



a— 161 — 1 



exp| it ^ XkUk I du\ ■ ■ ■ du r , 



expjzi J x(p{dx]^j, 



where (p is a Dirichlet process with parameter a defined by a{xk} = bk, k = 1, 
...,n, and a{0} = a — |b|. Hence, n $ can be viewed as the characteristic 
function of the mean of a Dirichlet process with a parameter a supported 
by a finite set. 
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The function n <3? admits an interesting representation, as a single contour 
integral, which Erdelyi achieved by resorting to the following argument. Take 
the Laplace transform of 



ip*(a) := (T(a)) a a n $(&i, b n ; a; itaxx, itax n ), 



namely, 



r+oo 

/ ip*(a)e- wa da 
Jo 

- w a r-+co f r -I 

/ a a exp< —wa + ita / xip(dx) > da 
X{a) Jo I J J 



(u; — it 



■ cxp 



log 1 



itx\ 

— ]aldx) 
w I 



[Re(w) > 0] 



[by Theorem l(iii) 



k=l 



itxk \ 
w J 



The well-known complex inversion formula for the Laplace transform [see, 
e.g., Henrici (1991), page 278] yields 



x 



a-l 



r(a) 



1 /-7+i-R 

lim / e 

i?^+oo 27TZ J^-iR 



11 1 

fe=l 



7>0, 



which, letting x = 1, gives 
^(e^W) - r(a) 



lim 

2iri R— >+ooj^_j^ 



e -i-n 1 



1(1" 



w 



dw 



PV / denoting principal value integral. Moreover, from the Cauchy integral 
theorem, the path of integration can be deformed into any contour which 
consists of a simple loop, (£t, beginning and ending at — oo, and encircling all 
the finite singularities of the integrand, that is: 0, itx±, . . . , itx n . See Figure 2. 

The above discussion leads to a first representation — under some restrictions- 
for the characteristic function of <p(I)- 
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Theorem 6. Let a be a measure with a bounded support S C R, and, 
for any t m R, let St be the closure of the convex hull of {itx :x £ S}. Then, 
with a = a(S) > 0, 

/ e ^(/)^ ( ^) = £Wp V r +lC °e-^e W \-c(--;a,l)\d W 

JP 2lTl V 7 -ioo ^ a I \ W J J 

/ioWs /or every 7 > 0. Moreover, if <&t is any contour consisting of a simple 
loop, beginning and ending at —00, and encircling St in such a way that 
S t n£ t = 0, then 

f e^(^) = ^ / ^exp{- C (--;a,/)W 
jp 2-iri Jt t w a I V w /J 

Proof. In view of the previous discussion, the statements are true if 
S is finite. More in general, if S is bounded, resort to the approximating 
functions £ n , used in the proof of Theorem 1, to write 

e u ^£> a (d<p)= lim / e ittp{ -^9Jd V ) 

n—i+oojf 

= V~^ / ^ ex p\-((- — '> a n,l)\ dw. 

ZTTl n-»+oo7 £t w a I \ W /J 

Next, interchanging limit and integral gives 

lij P / -r ex p(-cf-— ;an,^l ^ = / -^ ex p( _ cf-— ;a:,l)\ dw. 

n^+ooj tt W a I V W /J j£ t W a { \ W J J 

Finally, by the Cauchy integral theorem, one has 

— exp<^ -C ;a,J \ dw = PV — -expl -([ ;a,l) \ dw 



tt W u l \ W J ) J-y-ioo W I \ W 
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for every 7 > 0. □ 

To deal with general parameters a, assume that condition (ii) in Theo- 
rem 1 is satisfied and, for any k in N, define as 

a {k) (-) = 6- k (-)a(-oo,-k] +a(-n (-k,k)) + 5 k (-)a[k, +00). 

It can be proved that the distribution of / x(pk{dx), where (pk is the Dirichlet 
process with parameter , converges weakly to the mean of the Dirichlet 
process with parameter a, as k tends to +00. Hence, by the continuity 
theorem, 

(7.1) / t^f^9 aW {dip)^ / e u ^@ a {dtp), iSl,fc^+oo. 
From Theorem 6 it follows 

/ e M^ )W = E&PV r +i0 °^exp{-c(-^),/) }dw, 

(7.2) 

which proves 

Theorem 7. Let a be a measure in ¥ a> for some a > 0. Then 
(7-3) JF 

= hm MpvT^expl-cf-^M)^ 
fc^+oo 27^^ J-y-ioo w a I V w /J 

holds for every t in R and for any strictly positive 7. 

Remark 3. Principal value integral in the right-hand side of (7.3) re- 
duces to Lebesgue integral if a > 1. This is the case when, in the presence 
of exchangeable observations, one is interested in the Fourier transform of 
any conditional p.d. of </3(I), given those observations. 

As a first application of this result, consider the problem of characterizing 
symmetric p.d.'s in M a for some fixed a > 0. It is well known that the p.d. of 
<f(I), under @ a , is symmetric if a is symmetric. See, for example, Regazzini, 
Guglielmi and Di Nunno (2002). Here, in addition to a simple proof of this 
fact, we prove the converse. 

Theorem 8. Let a be an element in ¥ a . Then the distribution [i a of the 
mean of a Dirichlet process with parameter a is symmetric if and only if a 
is symmetric. 
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Proof. Suppose that a in F a is symmetric. Then, for any fcgN, 
r(a)pv ^'»^ exp |_ I ^ ^ 1 ^ 



2-iri 7 7 - loo ^ a I Vr V w J 
T(a) n^ 00 e w 



2iri 



'7—100 

which, by the change of variable y = —x, can be shown to be equal to 

2lTl Jj-ioo w a I Jr \ w J J 

thus implying, by virtue of Theorem 7, that the characteristic function of 
if (I) is real- valued. 

Conversely, if (i a in M a is symmetric, one has 



l im Mpv / 7+80 °£!lexp(- / logfl-— )d(a-A^(-x))\du 

^+00 2lll J-y-ioo W a { JWL \ W / J 

= l im iMpv r +l °°^-exp(- [ \oJl + —\dA { - k \x)\d W 

k->+oo 2-Kl Jy-ioo W a { Jr \ W J ) 



(by change of variable) 



lim £W PV r + - ^ r_ f _ ^\ dAik){x) 

k^+oo 2-Kl J-y-ioo W a I Jm. V W J 



dw 



(by symmetry of fi a ). 
The proof is concluded by resorting to Theorem 2 (uniqueness). □ 

Another application concerns the problem of determining the p.d. of 
</?(/) = f f d<p, when dp has the functional Dirichlet distribution 3) a , and 
/ is any measurable real-valued function. Here, Ql a can be understood as 
the functional Dirichlet distribution of a random probability measure fh 
on a probability space (f2,J^); that is, £1 and & need not coincide with 
M and ^(M), respectively. In this general framework, a must be a finite 
measure on (fi, J^) satisfying a = a(ft) > 0. If / is bounded, then the dis- 
cretization process used to prove Theorem 6 can be applied to give 

9a{e itrntf)) = £W PV / 7+l0 °^exp{- / logfl-^ao /"!((&)} dto, 

2lTl J 7 -joo W I JR V W J J 

allowing an obvious generalization of Theorem 8. 

Theorem 9. Let m be a random probability measure on (f2,J?) with 
distribution 3l a , and let f be a measurable function from £1 to ~R such 
that J^log(l + |/|) da < +00. Then, the p.d. of fh{f) is the same as the 
p.d. of (p{I), where p is a random probability measure on (R,^(R)) with 
p.d. S> a0 f-i. 
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8. Vector of means of a single Dirichlet process. Let (Q, J£",P) and Ql a 
be those defined in the final part of the previous section, and let fi,---,fd 
be measurable functions from Q to R satisfying 

(8.1) I log(l + |/ fc |) da < +oo, k = l,...,d. 

Jn 

For any t = (t±, . . . , td) in R rf , define a{t,f) [with f = (/i, . . . , fd)] as the im- 
age measure — defined on (R,^(R)) — of a, through the measurable function 
(t,f) : V'l i './//• that is, 

«(t,f> = "° (t,f) _1 - 

If m is a random probability measure on (Q, J£") with distribution 2l a , then 
Theorem 9 states that the p.d. of 

^tjfj dm 

coincides with the p.d. of <f(I), where (p is a random probability measure 
on (R,^(R)) having p.d. @ a{t f) . 

8.1. Multidimensional Lauricella identity. These remarks, combined with 
a straightforward application of Theorem 4, yields the following multidimen- 
sional form of the Lauricella identity, that is: Suppose (8.1) holds true, with 
a = a(£l) > 0, and denote the p.d. of (m(/i), . . . ,fh{fd)) by // Q) f . Then, for 
any t = (t\ , . . . , td) inW 1 , one has 

f 1 7 rT-A*af(^ x ) 

V (l + i(t,x)) cPQ,U ; 

= / exp< — / log(l + w(t,f)) da >B( du;c, a — c) 
J [o,i] I J 

if a> c> 0, and 

L (i+i(t,x))^- f(dx) 

_ r(c-a + l)r(a) 
~~ 2mT(c) 

x J^ +] expj- y log(l + itu(t,f)) daj™ c_1 (™ - l)"-^ 1 dw 

if c> a > 0. 

This proposition has also been proved in Regazzini, Guglielmi and Di Nunno 
(2002) when a = c and in Kerov and Tsilevich (1998) when a = c = 1. The 
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following sections illustrate two applications. The former concerns the proof 
of absolute continuity, with respect to the Lebesgue measure on R. d , of /x Qj f . 
Regarding this point, Firmani (2002) has determined a density function of 
H a j by inversion of the above identities with d = 2, and for suitable choices 
of fi and /2 . The latter application deals with the problem of determining an 
expression for the moment generating function of the variance of a Dirichlet 
process. 

8.2. Absolute continuity of fj, a f. In view of Proposition 2 and of the 
remarks at the beginning of this section, the p.d. of m((t,f)) is absolutely 
continuous with respect to the Lebesgue measure on K if f := (/i, . . . , f d ) 
is not affinely a -degenerate; that is to say, there are no v in M. d \ {0} and 
b in M for which a{(v,f) = b} = a. As far as absolute continuity of // Qi f is 
concerned, the following statement is valid. 

Theorem 10. Suppose a and f satisfy (8.1) and assume that f is not 
affinely a-degenerate. Then the p.d. fi a j of (m(/i), . . . , is absolutely 
continuous with respect to the Lebesgue measure on M. d . 

Proof. The contrapositive statement will be proved. Indeed, if \i a $ 
is not absolutely continuous, there are e > and, for every 5 > 0, a set 
A = UfcLi A fc with Ai = x d =1 {ai i] ,bi l) ), Ai n Aj = if i / j, and X d (A) = 

E¥=iU.t=i(.b ( r k) - a ( r k) ) < 5, such that a {m(f) £ A} = fi a ,f(A) > e. Denote 
the Lebesgue measure on M. d by X d and choose t\,. . . ,td in such a way that 

<tj < mm 1 < k < M Uie{i,...j-ij+i,...4}{ b k ~ 4 )M Then > 

M d M d 

k=ij=i k=ij=i 

and, therefore, letting A' = {J^ =1 faa^ + ■■■ + t d a ( k d) ,t 1 b£ ) + --- + t d b^ ] ) , one 

has A 1 (A') < J2kU Ej=i tj{b k j) ~ a?) < 5 - Moreover, A C {x G M d : (t, x) G 
^4'} and e > is such that, with each 5 > 0, one can associate a set A' G =5^(R) 
satisfying A 1 (yl / ) < 5 and £^ Q {m((t,f)) G A'} > e, contradicting the absolute 
continuity of the p.d. of m((t,f)). □ 

8.3. Moment generating function of the variance of if . If a satisfies con- 
dition (ii) in Theorem 1, then the random variance 

V:= [ {x-v{I)f(p{dx) 

is finite, &.s.S a . Some of the results given in Section 7 are employed here 
to determine an expression for the moment generating function, gy, of V . 
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The definition used for g v coincides with the one formulated, for instance, 
in Section 13.5 in Fristedt and Gray (1997), so that 

(8.2) g y {t) = @ a {e-*) = ® a (exp{-t(<p(f 2 ) - (<^(/i)) 2 )}), * > 0, 

where fi(x) = x and f2(x) = x 2 for every x in R. After observing that e tx2 
is the moment generating function, at x, of a Gaussian-distributed random 
variable with zero mean and variance equal to 2t, apply Fubini's theorem to 
obtain 

If the support of a is bounded, by using arguments similar to those employed 
in the proof of Theorem 6 and by resorting to Theorem 9, with / = ufi — tfy, 
it is straightforward to prove that 

^f e Ulp(fl)-tip(f2)\ 



2iri 

holds, provided that 7 W $ is any element of (0, +oo) satisfying 

J u ,t> sup (ux — tx 2 ). 

xGsupp(o) 

The resulting expression for the moment generating function of V is 



4^3/2, 



f e -2/ 2 /4p V r»- 4 l ™?l ex J_ c (l.^ tf2 _ y ^ tf M dzdy 

m Vv, -ioo z a I \£ J) 



'y,t 

and 7^ is any element of (0, +oo) satisfying 

r y'y t t> SU P (yVtx-tx 2 ). 

xSsupp(o) 

An extension of this representation to cases in which the support of a is 
arbitrary can be obtained by considering the sequence of truncated measures 
(&^)k>i, as defined in Section 7. Since the distribution of the variance of the 
Dirichlet process with parameter converges weakly to the distribution 
of the variance of the Dirichlet process with parameter a, as k tends to +oo, 
a continuity theorem for moment generating functions yields 

r(a) 



4-7T 3 / 2 ? 



x lim / e- y / 4 

k— »+oo , 



x PV jT v '^J ^ exp|-C Q; , </ 2 - yVt/i) } dzdy. 



See, for example, Fristedt and Gray [(1997), page 262]. 
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